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Abstract 

Secondary reactions induced by relativistic beams in inverse kinematics in a thick 
target are relevant in several fields of experimental physics and technology, like sec- 
ondary radioactive beams, production of exotic nuclei close to the proton drip line, 
and cross-section measurements for applications of spallation reactions for energy 
production and incineration of nuclear wastes. A general mathematical formulation 
is presented and successively applied as a tool to disentangle the primary reaction 
yields from the secondary production in the measurement of fission of a 238 U pro- 
jectile impinging on a proton target at the energy of 1 A GeV. 

PACS: 25.40.Sc; 25.70.Mn; 25.85.Ge; 25.85.-w; 24.10.-i; 29.25.-t; 29.30.Aj 
keywords: SECONDARY AND MULTIPLE REACTIONS IN INVERSE KINE- 
MATICS; In-flight identification in Z and A by magnetic spectrometer, Measured 
primary spallation and fragmentation cross sections; Fission and evaporation residues 
cross sections; Spallation reaction, p( 238 U,X), E=l AGeV; Relevance for the pro- 
duction of radioactive beams. 



1 Introduction 



The reaction residues of projectile spallation and fragmentation, produced 
with heavy-ion beams can undergo consecutive nuclear reactions in the tar- 
get. The secondary reaction production, especially by the use of a thick target, 
can extend towards exotic nuclei, which would not be generated in the primary 
reaction. The induction of multiple reactions in thick targets is also a technique 
to obtain exotic beams. On the other hand, in the measurements dedicated 
to the extraction of formation cross sections, the secondary reactions enter 
in competition with the primary production and are a disturbing process. In 
this work we will refer to a recent measurement of the isotopic production 
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cross section in the reaction 238 U+p at 1 A GeV [1,2,3]. The experiment was 
performed in inverse kinematics with the high-resolution magnetic spectrom- 
eter FRS [4] (GSI, Darmstadt) directing the uranium beam on a hydrogen 
target. The target thickness was optimised in order to maximize the primary 
production rate in respect to parasitic reactions and to limit the average share 
of secondary reactions to values of about 5%. This share and its uncertainty 
enters in the total experimental uncertainty; though, some of the measured 
isotopes, suffer from a much larger competition of the secondary reactions. 

When the primary isotopic production extends over a wide neutron-number 
range, as it will be pointed out, we can not extract reliable cross sections in 
the neutron-deficient side. Especially for the fission residues, the approaching 
towards the residue corridor defines a limit for the measurement technique: 
in this case the neutron-deficient isotopes will reveal to be mostly or entirely 
produced by secondary spallation of fission residues. 

The beginning of the present work is dedicated to the derivation of the mathe- 
matical formalism to treat consecutive reactions in a general form, indepen- 
dently of any application. Furthermore, we will focus on the cross-section mea- 
surement, and an analytical recipe to disentangle the secondary contribution 
from the primary production is presented. The method has been successfully 
applied in the data analysis of the reaction 238 U+p at 1 A GeV. 



2 Derivation of the exact formulation for multiple reactions 

2.1 Attenuation of the beam 

When a beam of ions characterised by atomic number Z and neutron number 
N interacts with a layer of matter, the intensity T evolves accordingly to the 
total reaction cross-section of the projectile o"o = o"(N Zo) and the properties 
of the target, like the density p[gcm~ 3 ] and the mass of the target nuclei 
Aar [<?]■ The loss — dX (£) of the beam intensity is proportional to the path d£ 
traversed inside of the target medium and to the intensity Xo(0 reached at 
the position £ in the beam direction: 

-dl (0 = a tar cr • 2o(Od£. 

The target properties are represented by a tar [cm~ 1 mb~~ 1 ] = 10~ 27 -J\f A ■ p/A taT , 
with M A indicating the Avogadro's number. Integrating over the whole path 
x we obtain the probability for the projectile (N Z ) to survive for a path x 
in the layer of matter as 

P (N Z ,x) = X (x)/T (0) = e-^tar* = e -a oXj (1) 
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where x [ m b *] — a t&rX indicates the number of target atoms per 10 27 cm 2 . 
2.2 Probability of primary reaction 

The probability for the projectile (NoZo) to interact with the target in a path 
length d£ with a total reaction cross-section a and produce a residue (NiZi) 
with a probability P(o->i) is defined as 

dP(N Z , NiZi, df) = atard^ • <?o ■ P(o->i) = d C ■ er(o-»i), 

where O"(o^i) = o"o • P(o->i) results to be the production cross section for the 
reaction (N Z — > NiZi), and d( = a tar d£ is the number of atoms in a layer 
of matter defined by the path d£ per 10~ 27 cm 2 . The total probability for the 
projectile (NoZo) to react only once in the target in any position £ and produce 
the observed fragment (NiZi) is expressed by the probability for the projectile 
to survive for a path £, react in d£ and produce a residue that will traverse 
the remaining length of the target without any further reaction; as pictured 
in the drowing below, this expression should be integrated over any path £ in 
the form : 

f [Po(N Z , • dP(N Z , NxZx, d£) • = 
Jo 

= a^e-^X J* p (N Z , C) e°< d(. 



N Z 




t x 

Introducing Eq. (1) in the integral, the solution expressed in terms of target 
thickness gives: 

/ e -croX e ~ ai X \ 
Pi (N Z , NiZi, X ) = -<T(o-+i) + • (2) 

V <7 - <7i G\-Oq) 

2.3 Secondary reactions 

After the primary reaction of the beam occurred, we should consider the pos- 
sibility for a further interaction of the residues with the target. In this case, in 
order to obtain the probability of observing a secondary fragment N 2 Z 2 , the 
quantity to integrate over any path through the target in beam direction is 
the product of three terms: the probability for the beam to undergo a primary 
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reaction during the path £ and produce the intermediate fragment N1Z1, the 
probability for a further reaction (N1Z1 — > N 2 Z 2 ) in the interval d£, and the 
probability that the final residue survives for the remaining path x — £. This 
description is represented below in the drawing and it results into the following 
equation: 



f 

Jo 



7MN Zo, N i z i> • dP(N 1 Z 1 , N 2 Z 2 , d£) ■ P (N 2 Z 2 , x-£) 
= <J(i->2) e~ a ^ j* Pi(N Z , NxZx, C) dC 



N n Z 



0^0 



NiZj 



T 



N 2 Z 



2^2 



X 



Introducing Eq. (2) in the integral we obtain for the secondary reactions: 



PaCNoZo.NxZ^NaZa,*) 



(7(0-^1)^(1^2) 



e -°oX 



+ 



e -ViX 



+ 



e -(?2X 



(<J -(7i)((Jo-(72) (<7i-<7 )(<7i-<7 2 ) ((7 2 -(Jo) ((7 2 -<7i) 



(3) 



As it will be demonstrated in the section A.l, the procedure followed so far 
to obtain the probability for secondary reactions can be extended to higher 
orders according to the recursive relation: 



P n (N Z ,NiZi,... ,N n Z n ,x) 



<y(n-l^n)e~ anX £ 7Vl(N Z , NiZi, • • • , ^n-lln-1, C) ^ d(. 



(4) 



Repeating the iterative integral for higher and higher orders we obtain the 
solution for the n th order: 



n 

^(NoZcNxZx,..- ,N n Z„, X ) = (-1)" n^-i^) E ~ 



i=0 I! (&i - 



(5) 



3 A more stable approximated formulation 



If we neglect the technical and chemical properties of the target, i.e. we con- 
sider the target homogeneous and a full acceptance of the spectrometer, the 
relation (5) is formally rigorous and exact. Nevertheless, it should be observed 
that the term (<7j — (Jj) could generate instability in numerical calculations: 
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when very few mass is removed in one step of the chain of consecutive re- 
actions, the total cross-sections <7j and Oj become almost identical, and the 
denominator (cr, — Oj) becomes very small. In this case, the sum in Eq. (5) 
remains finite but its numerical computation can be problematic in this form. 

A possibility to remove this inconvenience is to perform a series of consecutive 
l st -order approximations in iterating the integration (4). We start applying 
the approximation to the Eq. (2), which describes the probability of primary 
reaction: 

e -o"oX _ e ~°iX 
Pi(N Z ,NiZi,x) = -o-(o^i) ■ 

(Jo - <7i 

Expanding the exponential to the 2 nd order in (tfo— we obtain the simpler 
relation: 

<Jp + <J\ 

TMNoZcNiZi,*) « <7 M xe 2 X . (6) 
Integrating Eq. (6) by the relation (4) we obtain the probability of observing a 
secondary product generated by a given intermediate fragment in the following 
form: 



C(0-vl)C(l-v2) e 



OiXl\-[— 

Jo 



Expanding to the 1 st order in | (a + <?i) x ~ <y 2X-, integrating, and isolating 
X 2 /2, we obtain the approximated form of the Eq. (3) in the 1 st order : 



2 <Jq + Q\ + cr 2 
P 2 (NoZ ,NiZi,N 2 Z2,x) ~ o-(o->i)0-(i_>2) ye 3 

Physically, this relation can be intuitively pictured as "one-third-target ap- 
proximation", as the argument of the exponential represents the attenuation 
of the beam, the primary and the secondary residue, respectively, when they 
cross a third of the target. 

Iterating the integration (4) followed by the first order expansion in respect to 
the argument of the exponential, we obtain the approximation for n consecu- 
tive reactions, equivalent to dividing the target in n portions of equal thickness, 
each one originating a successive reaction with a probability a^-i^x/i' 



En _ 

P n (N Z , NxZi, • • • , N n Z„, x ) = f[ ( a(M) X ) e n+1 X 

i=l V % J 

y — — X 

= n^-i^)\e n + 1 . ( 8 ) 
i=i n! 

To obtain the Eq. (8) we can follow the prescription presented in the section 
A.2. 
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4 Measured reaction probability 



Besides the production rate in the target expressed in Eq. (5) , the measured 
production rate of an isotope NZ should also take into account the probability 
that the residue NZ is observed outside of the target, and is detected in the 
spectrometer. The measured probability to produce NZ can be written for a 
given chain of intermediate products as 

M n {N Z , • • • , N n Z n , x) = ^n(N Z , • • • , N n Z n , X ) ■ *n(NiZ 1? • • • , N n Z n ), 

where t n is the transmission coefficient depending on the kinematics of all the 
intermediate isotopes: it depends on the instrumental device and it represents 
the probability that a fragment with a given velocity is transmitted through 
the spectrometer. The transmission is a key parameter in the measurement 
of fission products which are only partially accepted, as some of the residues 
are emitted with large angle and hit the pipe of the spectrometer. The mea- 
sured probability to obtain NZ by multiple reactions accounts for the sum of 
the orders of the reactions, that is the superposition of the contributions of 
primary, secondary and multiple reactions evaluated for any possible choice of 
intermediate fragments 

M(N Z ,NZ,x) = E E A4n(NoZ ,---,N n Z n ,x), 

order n n—l chain 

where the first sum accounts for the order of reactions and the second sum 
accounts for an order given for all the possible chains of successive n — l 
intermediate fragments leading from the initial projectile N Z to the observed 
residue NZ. We can express the whole relation in terms of cross-sections if we 
introduce the apparent cross section 

1 A4(N Zq,NZ,x) 

^(N Z„- > NZ)(X) - , 

where t(NZ) is the transmission factor of the spectrometer. It is derived from 
the measured spectrum of longitudinal velocities, by assuming the emission 
isotropic in respect to the centre of mass; this assumption allows to evaluate 
the angular distribution of the fragments and the fraction which is selected 
by the angular acceptance. Due to the difference in the distributions of the 
emission velocities related to fission fragments or evaporation residues, re- 
spectively , the coefficient t(NZ) depends strongly on the dominant reaction 
process responsible for the production of a given isotope NZ. In order to have 
a more complete description of the reaction process, we can disentangle fission 
residues from evaporation residues. If we assume that fission could occur only 



6 



once in a chain of successive reactions we obtain the general description: 
1 



°"(N Z ->NZ)(X) - - E E 

A- order n n—1 chain 



with: 



1Z cvr = 



i=NZ 



n «ux) 



=NiZi 



j=NZ / evr, evr, 
fisj SJ- fisj 
n 



j=NZ 



'Ml 



(9) 



^ = ~y jX II (°ZUx)> 



j=NiZi 



/T-cvr 
1 n 



C r (N 1 Z 1 ,N 2 Z 2 ,--- ^^xZ^x.NZ) 



t(NZ) 



e fi v s r ; C SJ (N 1 Z 1 ,N 2 Z 2 ,-- - .N^tZ^.NZ) 

t(NZ) 



^(NoZo.NiZi,... ,N n Z n ,x) = (-1) 



71! 



E 



n + 1 



o n [pix - o-jx) 



The first sum describes the apparent cross-section as a composition of con- 
secutive orders of multiple reactions. The second sum over all the possible 
intermediate fragments can be represented by an expansion in (n — 1) chains, 
the i th one satisfying the condition (Nj_i + Zj_i) > (Nj + Zj) > (N i+1 + Z i+1 ): 
each successive intermediate fragment should in fact have lower mass but not 
necessarily lower charge due to the charge exchange. Nevertheless, since this 
process has very low cross section, the condition can be safely decomposed 
into Nj_! > Nj > N i+ i and Z*_i > Z { > Z i+1 for each sum i. K e J r and 7^ vr ' fisi 
are the terms defining the reaction chain: is the term representing the se- 
ries of reactions producing NZ as a chain of consecutive evaporation residues; 
7^evr, fisj differs f rom the former one because at the step j a fission reaction 
occours. Due to the very low fissility of the fission residues, the probability to 
undergo two consecutive fission reactions is negligible: therefore, we imposed 
only one possible fission event in the chain. The terms 7^ evr and 7^ cvr ' fis i are 
transmission coefficients ratios and depend on the optics of the instrumental 
device; a more detailed description of these coefficients will follow in the next 
paragraph. A n is the term describing the attenuation of each nucleus travers- 
ing the target, either the beam-like projectile and the consecutive residues. It 
should be observed that this terms takes into consideration also the attenua- 
tion of the final residue, due to the probability to undergo a reaction of order 
7i+l; nevertheless no reaction of order n + 1 appears in the term describing 
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the reaction chain. As discussed in the previous paragraph, according to the 
possibility to have an unstable solution, it is important to decide in which 
condition the exact form containing the term ai — aj and referred to in Eq. (5) 
can be applied; it could be safer to use the approximated form represented by 
the Eq. (8) which remains valid as long as the products (crj — <Tj)x are small 
compared to unity. 

We simulated the secondary and ternary production (expressed as apparent 
cross-sections) in the reaction of 238 U at 1A GeV impinging on a target of hy- 
drogen of 87.3 10 -3 [g/cm 2 ] by solving the system (9) with the choice of the ap- 
proximated form (8). The production cross-sections c^I Zo -*n-z- an d ct n z -»nz- 




Fig. 1. Top. Experimental cross-sections for the reaction U+p at 1 A GeV 
[1,2,3] measured in inverse kinematics and corrected for the secondary contribution. 
A hydrogen target of thickness 87.3 rag j err? was used. Bottom. Contribution of 
the secondary reactions to the appearent cross section as calculated with Eq. (9). 
The secondary production is more than one order of magnitude lower than the 
primary production. The secondary fission residues, shown on the left, account for 
primary evaporation residues followed by secondary fission residues or primary fis- 
sion residues followed by secondary evaporation residues. The primary evaporation 
residues followed by secondary evaporation residues are shown on the right. 
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were measured experimentally in inverse kinematics at the FRagment Sepa- 
rator in Darmstadt [1,2,3]. The total cross-sections <7j have been calculated 
numerically according to the model of Karol-Brohm [5,6], very well suited for 
proton reactions. Another choice would be the formula of Benesh et al. [7], 
where the total cross-section is calculated neglecting the energy of the projec- 
tile and the nuclear properties of the target: as a consequence of these approx- 
imations the description of ion-proton reactions gives an almost-systematic 
over-prediction of around 20% in the total cross section; nevertheless, in the 
case of ion-ion reactions, the formula of Benesh et al. gives satisfactory results 
and is then preferable due to the lower computing time. The fundamental 
parameters to enter in the system of Eq. (9) in order to obtain a physical 




Fig. 2. Study of the evolution of the mass distribution of the residues generated 
in the reaction 238 U+p at 1 A GeV for different orders of multiple reactions in a 
hydrogen target of thickness 87. 3mg/ cm 2 . The primary production corresponds to 
the experimental (dots) cross sections [1,2,3] corrected for the secondary contribu- 
tion. Nine elements from Tb to Ta are still not analysed and their cross sections are 
replaced by a calculation performed with the codes BURST [8] and ABLA-PROFI 
[9,10,11] coupled together. The secondary and ternary production is calculated ac- 
cording to Eq. (9). The calculation has been performed using both the exact form 
(dots) of the term A n and the approximation of equally-partitioned target (his- 
togram) described in Eq. (8); the perfect overlapping of the two calculations assures 
the reliability of the approximation. 
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description of multiple reactions are the production cross sections of"^ and 
offj^j. Due to the enormous number of these parameters, a very fast numerical 
calculation is needed. The best solution we found is to couple a parameterisa- 
tion to obtain the hot prefragments with a physical fission-evaporation code. 
The routine BURST (applied for the correction for secondary reaction in the 
analysis presented in [8]) was chosen to reproduce the intra-nuclear cascade. 
The evaporation-stage has been simulated with the statistical de-excitation 
code ABLA [9,10] and the fission code PROFI [11]. Also the possibility of 
multifragment emission, in the case of light and highly excited systems, has 
been taken into account and simulated according to the reference [12]. The 




Fig. 3. Evolution of the isotopic distributions for different orders of multiple re- 
actions in the reaction 238 U+p at 1 A GeV. The different productions have been 
obtained as discussed in fig. (1) and in fig. (2). Evaporation residues (top): The 
mean neutron number of elements close to the projectile like Fr are not affected 
by consecutive reactions, while a shift towards lower neutron-numbers appears for 
elements like Os, produced in the light side of the mass-distribution slope. Fission 
residues (bottom): The isotopic distribution of multiple-reaction products is de- 
formed and the mean neutron number tends to move towards the residue corridor 
(dashed line). 
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resulting isotopic production probability has then been normalized to the total 
reaction cross section obtained with the model of Karol-Brohm. 

We verified that the results fit well with the available data of proton-reactions 
for 238 U [1,2,3] and 208 Pb [8] at 1 A GeV and 197 Au [13,14] at 800 A MeV. 
In fig. (1) the result of the calculation for 238 U+p at 1 A GeV is shown. The 
primary production is represented by the experimental cross sections; due 
to the competition between fission and evaporation residues, nine elements 
(from Tb to Ta) have not yet been analysed; for the calculation, this gap 
has been filled with model-deduced cross sections. It should be noted that in 
this figure we anticipate the result on the experimental primary production 
cross section, that have been disentangled from the secondary contribution 
by applying the procedure described in the following section. The secondary 
yields are one to two orders of magnitude lower than the primary production as 
shown in fig. (2). The secondary evaporation residues are expected to extend 
to lighter masses and slightly overlap with the heavier fission products. In 
general, even if the mean neutron number of the isotopic distribution does 
not change significantly in respect to the primary evaporation residues, we 
observe an increased production of neutron-loss channels and an extension 
of the isotopic distribution towards the proton-rich side. In contrast to the 
case of evaporation residues, the secondary fission residues are evidently less 
neutron rich. This is almost exclusively due to secondary evaporation residues 
that reduce the primary heavy neutron-rich fission products and populate the 
neutron deficient side; this is shown in fig. (3). The ternary production, three 
to four orders of magnitude lower than the primary production, shows an 
even more enhanced production around the residue corridor. The production 
resulting from secondary and ternary reactions, could be increased using a 
thicker target than the one considered in the present calculation. 



5 A recipe to disentangle primary and secondary reactions 

The calculation presented in fig. (1) derives from the direct solution of the 
system of equations (9), applied to the measured reaction cross section of 
238rj_|_p a t 1 A GeV. Evidently, the result of the experiment was the cumulative 
detection of the primary production together with the secondary, without 
any clear experimental indication about how to select the primary yields. 
Thus, the primary-reaction cross section were disentangled from the secondary 
contribution solving the system (9) inversely. 

From the calculation shown in fig. (2), we confirm that the ternary contribu- 
tion induced by a thin target is even lower than the uncertainties of the data 
introduced by the measurement, and it can be safely neglected. In this case, 
the notation has been changed in respect to Eq. (9). The index "fis,evr" defines 
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primary fission residues leading to secondary evaporation residues, "evr,fis" de- 
fines primary evaporation residues undergoing secondary fission and "evr,evr" 
primary evaporation residues producing secondary evaporation residues. This 
consideration justifies the reduction of the system (9) to the second order of 
consecutive reactions only. Thus, we can write the matrix relation 



^£>fis /2^ s X^ S -\- Q^ s ' evr ^ evr _|_ (^cvr,fis^>fis 



(10) 

(11) 



where X is the vector of apparent cross sections directly measured by the 
experiment and X is the set of the unknown primary cross sections that we 
intend to extract: 



X 



The elements 



C(NoZo^NoZo) (x) > 
C(N Zo-No-iZ ) (x) 

°"(N Z ^NZ) (x) 



V 



X = 



1 C(N Zo^N Zo) 
C(N Zo->No-iZo) 



C(N Zo^NZ) 



«nz = - ^i(N Z ,NZ,x) 
Pn^ = \ a (NiZi _ >NZ) x^ 2 (N Z ,N,Z J ,NZ,x)r 2 (N,Z l ,NZ) (12) 

are collected into the diagonal matrix C and the triangular matrix Q respec- 
tively : 



a N n _,Zn 



/ /3n z 
/?n z /3n _iZ 



, Q 



\ 



/3n z 



/3n,z+i • • • /3n+i,z /?nz 



V 



/ 



The terms /^Z; are evaluated for any of the involved secondary reaction 
through a model calculation of each intermediate cross section (X^z^nz)- The 
attenuation terms A\ and A2, contained in a^z and in ftn^, respectively, are 
obtained by the calculation of the total reaction cross sections cXj through a 
model. The transmission coefficient ratios 7^(NjZj,NZ) are calculated for any 
combination of intermediate residue NjZj and final residue NZ as explained at 
the end of the present section. 
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The Eq. (10) describes the measured evaporation-residue yields as the sum 
of two contributions: the quantity £ evr represents the loss of primary cross- 
section due to the attenuation in the target. The term g cvr > evr is the gain factor, 
which takes into account the increasing of the yields due to the residues pro- 
duced in two consecutive evaporation-residue steps. The Eq. (11) refers to the 
fission-like yields, i.e. the production-rate of those residues which the measure- 
ment attributes to fission events according to the kinematical identification. 
The apparent fission cross-section suffers from the attenuation in the target, 
expressed by the term £ fis , and from two sources of gain: secondary evapora- 
tion residues of fission residues, represented by CJ fis ' cvr and evaporation residues 
followed by fission events, as described by CJ cvr ' fis . In the case of evaporation 
residues, the loss term depletes the yields of the elements close to the projec- 
tile. Part of these losses are redistributed by the gain factor Q evr and populate 
the tail of light evaporation residues. The effect of this process results into re- 
ducing the slope of the mass distribution of the evaporation residues. Another 
part of the losses goes into the term C7 evr > fis ; establishing a coupling between 
the Eq. (10) and the Eq. (11) and populates the fission yields with secondary 
fission fragments. Since the primary evaporation residues that could undergo 
a secondary fission event should be rather close to the projectile, the reac- 
tion products should preserve essentially the structure of the primary fission 
residues and do not introduce any modification or displacement in the origi- 
nal primary distribution. This is no more true when the primary reaction is a 
fission event. In this case, the losses are redistributed by the gain factor £ fis > evr 
as secondary evaporation residues of primary fission residues. The effect of 
this process is a strong deformation of the fission-fragment distribution. The 
primary fission products which, in general, are neutron-rich are turned by sec- 
ondary reactions into evaporation residues, generally neutron deficient, and 
have the tendency to cumulate around the residue corridor. As a consequence, 
the whole fission distribution appears less neutron-rich in the measurement 
than it is expected in reality. Moreover, the steeper is the slope of the fission 
distribution in the neutron-deficient side, the higher is the probability that 
some residues are produced by secondary reactions in that region. 



An example of the reaction mechanism ending up into the production of a sec- 
ondary fission residue is illustrated in fig. (4), where the isotopic distribution 
of the formation cross section of 110 Pd and 100 Pd for each possible intermedi- 
ate mother nucleus is presented. We deduce that the stable and neutron-rich 
110 Pd is produced with high cross section mostly by primary fission of heavy 
fragments, and its formation as a secondary residue should not enter in com- 
petition with its production by primary reactions. Inversely, the secondary 
formation of 100 Pd, close to the residue corridor, results to depend almost 
exclusively on secondary evaporation of the primary fission fragments: the 
measured yield of formation of 100 Pd should then be rather different from its 
primary production cross-section. 
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Fig. 4. Representation of the intermediate fragments NjZj that, once produced as residues of the reaction U+p at 1 A GeV, 
react a second time in the target and produce either 100 Pd (left) or 110 Pd (right). On the isotopic chart we represent the secondary 
production cross sections for all the possible reactions NjZj — > 100 Pd or NjZ, — > 110 Pd, multiplied by ct238tj_,,nz- an d divided by the 
total production cross section for the primary residues of 238 U. As clearly shown in the Z-projection, while the secondary production 
of neutron rich 110 Pd can be equally attributed to fission or to evaporation residues of the neighbouring elements, the secondary 
production of the less neutron rich 100 Pd can be prevalently attributed to evaporation residues of elements lighter than Nd. 



Thanks to the triangularity of the matrix Q, the system can be solved itera- 
tively by the recurrence relation: 



r (N Z ^NZ) 



evr 



~evr \ ^ ^-evr oevr 

a (N Zo^NZ) ^(NoZo-vNfZijPNiZi 



(13) 



A similar equation holds for ct^n z ^nz)- Introducing the expressions (12) in 
Eq. (13) and substituting the term A2 in the approximated form we obtain 
the recursive relation: 



r (N Z ^NZ) 

-- E 
2 z - 



a (N Z ^NZ) 



" N,->N 



+ a 



-evr 
= CT, 



(N ( ,Z -NZ) 



TTl^NoZo + 0"Nz) , 



T evr fM 7 \ "/ t „fi ^ N ° Z ° ~ Z + ^ 

<J (N Zo^N l Zi) <T (N i Zi- > NZ) i 2 l^i^ij ^ 



n 4 >n ,_ 
z,>z 



r (N„Z ^NZ) 



^iNoZo^NiZO^NiZi^NZ)^ " (NjZj, NZ) + 



-(cnoZo + <7nz) 

+ 

fis, 



(NoZo^NiZ^^NiZ^NZ)^ S (N;Zj, NZ) 



— (°"NoZ — Z&NiZt + C"Nz) 

eo 



(14) 



The relation (14) can be solved numerically for decreasing masses of the ob- 
served fragments NZ. Following this order, the unknown primary reaction cross 
that appears in the system has been already calculated in 



section a 



evr, fis 

(NoZo^NiZi) 



the previous step of the iteration. 

The solution of the system (14) gives directly the value of the primary produc- 
tion cross-sections. As discussed in the references [8,14,15], a clear indication 
of the incidence of the secondary reactions on the experimental measurement 
can be deduced studying the variation of the correction factor 

^(NZ) = O"(N Z -^NZ) / 6"(N Z ->NZ) 

that we should apply to the apparent yields in order to obtain the primary 
cross-sections. In fig. (5) we mapped the values of /i(nz) for each measured 
isotope. The distribution shows mainly a plateau for a positive correction of 
around 5% to 7%: these values, originating from the first exponential term in 
Eq. (14) should compensate the effect of the losses. The homogeneity of the 
correction is then perturbed by the effect of the gain factor. In some cases, 
like the evaporation-residue production, this effect does not present complex 
features: generally, the secondary evaporation residues with low mass-loss in 
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Fig. 5. Bottom. Representation of the correction factor fo(NZ) f° r secondary reaction 
on an isotopic chart for each residue measured in the reaction 238 U+p at 1 A GeV 
[1,2,3]. Top. Isotopic distributions of /i(nz) f° r the elements Ge, Zr, Pd and Te, pro- 
duced by fission and Pb and Fr, produced as evaporation residues. The dashed lines 
represent the residue corridor, calculated with the physical reaction code ABLA, 
imposing the sequential evaporation of higly excited heavy prefragments. The cor- 
rection factor drops to very low values around the residue corridor: the steep slope 
of the "hook" shaped /i(nz) distribution determines the limit for the observation 
of primary reaction products and for the measurement of their isotopic cross sec- 
tions. The correction for 100 Pd and 110 Pd, indicated by arrows, can be compared to 
%• (4). 



respect to the projectile should distribute on the same ridge of the primary 
production without showing a sensible deviation; as a consequence, the iso- 
baric distributions should simply be scaled in height according to a constant 
factor and preserve the shape unchanged. This prescription could be followed 
as a rather good approximation in the case of heavy evaporation residues in 
reactions where no isotopes with high fissility are produced: the correction 
ft-(NZ) would reduce to a function of the mass-loss only. This was the approach 
applied in the analysis of the evaporation-residue production in the reaction 
of 197 Au+p at 800 MeV [13], for the reaction 208 Pb+p at 1 GeV [8], and for 
the reaction 208 Pb+d at 1 GeV [15]. For these systems, and for the analysis of 
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evaporation residues, this method would even have an advantage in respect to 
the application of the recursion (14): in fact, the assumption of an exclusive 
dependence of the correction fo(Nz) on the mass-loss and the total cross-section 
of the intermediate fragment cr( NjZj ) allows to build a parameterisation based 
on the experimental data. On the other hand, the relation (14) suffers from 
the uncertainties introduced by the theoretical calculation of the intermediate- 
fragment cross-sections o"(N;Z;^nz)- Indeed, for systems like 238 U, dominated by 
a complicated competition between the fission and evaporation-residue chan- 
nels, the correction factor /i(nz) can not be easily related to the mass-loss. 
Moreover, the distribution of the evaporation residues of 238 U+p widens con- 
siderably around Pb: as a consequence the rain of secondary residues can 
not cover homogeneously the whole isobaric distribution and it will populate 
prevalently the neutron deficient side. As shown in fig. (5) this effect appears 
in the neutron-distribution of the correction factor for the element Pb and it 
induces a slight reduction of fo(NZ) for low neutron numbers. This effect tends 
to disappear for heavier elements like Fr, characterized by an almost constant 
correction factor. 

The use of a detailed description of the secondary reaction mechanisms be- 
comes unavoidable when fission fragments should be treated: in this case fo(NZ) 
largely chang function of the neutron number and the relation (14) can 

not be substituted by a parameterisation. As firstly formulated and demon- 
strated in the analysis of the fission residues in the reactions 208 Pb+p at 1 
GeV [8] , the correction factor is characterised by a "hook" shape and reduces 
steeply when it approaches the residues corridor. However, the steepness of 
the "hook" reflects also the steepness of the neutron-deficient side of the yield 
distribution. In fig. (5) an increased isotopic dependence of /i(nz) is shown for 
decreasing element numbers. The isotope 110 Pd, whose formation mechanism 
has been illustrated above, is on the plateau of the /i(nz) function and the re- 
lated cross-section should be corrected almost exclusively for the losses; on the 
contrary, 100 Pd is located on the slope of the correction factor. Other elements 
like Ge or Zr show corrections that even approach 100% for their lightest 
isotopes: of course, when the correction is too large, the rejection of these iso- 
topes becomes necessary. Isotopes showing correction factors smaller than 0.5 
are omitted in the tabulation of cross sections in ref. [1]. The steep slope of 
the correction factor becomes a technical limit that excludes the possibility to 
extend the primary fission cross-section measurement to the neutron-deficient 
side. 

We should still point out that, as shown in the definition (12), the gain factors 
include still another term: the transmission ratio for the secondary fragments 
72(NjZj,NZ) should be taken into account. As shown in fig. (6), the omission 
of the transmission ratio leads to an overestimation of the correction factor. 
The geometry of the spectrometer, as discussed in [16], has to be evaluated 
precisely in order to determine transmission values, which are deduced relating 
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the experimentally observed velocities of the fragments, the acceptance of the 
spectrometer, and the beam velocity. The details of the calculation of these 
coefficients, used in this work, is described in [1]. In the case of intermediate 
fragments NjZ,, that are not detected, the transmission factor can not be 
measured directly. However, in this case we can still estimate the transmission 
factor because the velocities are known. In order to reduce the transmission 
ratios to only experimental values it is worth while applying the following 
approximations: 

TT^N^NZ) = ^'^(NZ 2 )^^ * X ' (15) 

rr fis (na,nz) = * r '*^ ,NZ) - 1, (16) 

T fi — fN-Z- NZ) - ^ ( NA ' NZ ) « ^ NA ) > 1 (17) 
T 2 (N.Z..NZ) - £(NZ) ~ tfis(NZ) >1, (17) 

The relation (15) follows from simple considerations on the angular distri- 
bution of evaporation residues. In case of double reactions, the variance of 
the angular distribution in the laboratory is simply the sum of individual 
variances, which in turn are proportional to each individual mass loss, im- 
plying that the total variance is proportional to the total mass loss either 
for a secondary or primary reaction. This leads to ^^(NiZi, NZ) » t^NZ) . 
when the primary production of NZ is mainly related to evaporation residues. 
We should observe that in some cases, when the primary production of NZ 
is mainly related to fission, and there is an additional contribution coming 
from secondary evaporation residues of primary evaporation residues, the ra- 
tio of Eq. (15) becomes greater than the unit. When fission is involved in 
the production of intermediate or secondary fragments, even though the same 




Fig. 6. Study of the contribution of transmission in the calculation of the correction 
factor hmz)- The transmission ratios are introduced as prescribed in Eq. (15,16,17) 
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addition rule for the variance is applied, the variance for a process leading 
to evaporation residues is assumed to be much lower than in the case of fis- 
sion, and it can be neglected. This leads to the relations (16) and (17). 
It should be observed that only the relation (17) is relevant in the calcula- 
tion of the correction factor /i(nz)- This derives from the observation that the 
variance of the angular distribution of a fission fragment NjZj is not strongly 
modified when the residue reacts again and produces a secondary residue by 
evaporation. We can then neglect the contribution of the secondary reaction 
to the variance and substitute the transmission coefficient t^' 6 "; that can not 
be calculated on the base of experimental observables, with the transmission 
coefficient t fis of a fission fragment NjZj; this leads to the the approximation 
t 2 s ' evr (NjZj, NZ) m t fis (NjZj) . Since the variance of the angular distribution for 
the fission fragment NZ is larger in respect to the lighter intermediate fission 
fragment NjZj, the resulting term "Zf s ' evr (NjZj, NZ) is larger than the unit. 



6 Conclusion 



This work presents a study on the formalism of multiple nuclear reactions 
induced in a target. A general model has been developed in order to simulate 
the isotopic yields of secondary reactions, distinguishing between fission and 
evaporation residues. The model described in this work became the base of the 
numerical code "SECONDARY" . The program has been tested to reproduce 
the secondary-reaction correction used in the analysis of the reactions 208 Pb+p 
at 1 A GeV [8] .The comparison resulted in good agreement with the previous 
calculations. The code SECONDARY, coupled with the reaction code ABLA 
[9,10] and PROFI [11] has then been applied in the final step of the data- 
analysis presented in [1], aimed to determine the isotopic cross section of the 
reaction 238 U+p at 1 A GeV. 
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A Mathematical appendix 



A.l Extension of the description of the secondary reactions to the n order. 



The key relation to extend the description of the secondary reactions to the 
n th order is condensed in the following equality: 



n— 1 i 

y — - — 

i=o n fa - o-j) 

3=0 

3±i 



1 

n-1 

n (<r„ - o-j) 

3=0 



(A.l) 



Eq. (A.l) is a particular case, obtained imposing k = n — 1 and A = a n , of the 
following more general form: 



k l 

<=0 (<Ti - A) I! fa - <Tj) 



3=0 



n(A-^) 



(A.2) 



3=0 



For fc = the equality is trivial, and it can also be checked for k — 1. To prove 
Eq. (A.2) for any order, we should verify that the expression is recursive. If 
Eq. (A.2) is true at the order k — 1, we can extend it to the order k with the 
following passages: 



£ i 

i=o (a, - A) n to - o-j 



jfc-1 

+ E 



CTfc — A 



afc A L n k - o-j) i=o - a) ft - ^) 



3=0 



3=0 
3& 



a k - A 
1 



fc-i 



fc-i 

n (<t*; - 

3=0 

37^ 



+ E 



- A 



fc-i 

+ E 



i 



3=0 



k-1 

E 



i 



;T ' X L U (cr fc - a,) *=o (^ - a fc ) ft 1 (a, - a,-) i=o (a, - A) n -<Tj). 

3=0 3=0 3=0 

37^ 3^i 3^< 
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or, applying Eq. A. 2 to the last two terms in the brackets: 



1 

a k — X 



1 



fc-i 

n (<7k - Oj) 

3=0 



+ 



1 



k-1 ' k-1 

n (<T»i - o-j) n (a - o-j) 

3=0 j=o 



1 



(A.3) 



n (a - a,) 

3=0 
3& 



The probability for n recursive reactions in the target, expressed in the Eq. (5) 
reduces to the Eq. (1,2,3) for n equal to 0,1, and 2, respectivly. We can demon- 
strate that it is generally true if we obtain it recursively from "P n _i(NoZo, NiZi, 
• • • , N n _iZ n _i, x) by applying Eq. (4). We reduce Eq. (5) to the order n — 1 
and we introduce it in Eq. (4): 



P^NoZcNxZx,... ,N„Z n , X ) = 



n—1 



i=0 



JO tt 



0"n - <7i)C 



3=0 
3^i 



i=o n (<Tt - 

j=o 



(A.4) 



If we apply the equality A.l we can prove that the Eq. (A.4) reduces to the 
form ( 5). 



A. 2 Extension of the approximated secondary reaction formalism to the n th 
order. 



To prove that the Eq. (8) is a recursive extension of the Eq. (7) we should 
demonstrate that it could be derived from the order n — 1, i.e. we have to solve 
the integral 



nc(i-i-n) e 

i — 1 



-o-nx r (eg 

Jo (n-1)! 



n 



-o- n )C 
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If, before integrating, we expand the exponential to the 1 st order in ^ Y^]=o a j(~ 
a n (, we reduce to the approximated form: 



ii^i-i-i) e 



-<?nX 



^n— 1 

n 



— a„ 



X 



n+l 



j=0 gj 



(n- l)!(n+ 1) 



" Y 

= n^-i^)^ e_<7nXe 



— 0"r7 



n + l 



n! 



where the resulting expression is equal to Eq. (8). 
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